Orbital ordering, Jahn- Teller distortion, and resonant x-ray scattering in KCuF 3 
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The orbital, lattice, and spin ordering phenomena in KCUF3 are investigated by means of LDA+U 
calculations, based on ab initio pseudopotentials. We examine the Cu-3d orbital ordering and the 
associated Jahn- Teller distortion in several different spin-ordered structures of KCUF3. The ground 
state is correctly predicted to be an A-type antiferromagnetic structure, and the calculated Jahn- 
Teller distortion agrees also well with experiment. Concerning the orbital ordering, we find that even 
for a highly ionic compound such as KCUF3, the orbital-order parameter is significantly reduced with 
respect to its nominal value due to Cu(3d)-F(2p) hybridization. We also calculate the Cu K-edge 
resonant x-ray scattering spectra for Bragg reflections associated with orbital order. Consistent with 
previous studies, we find that the resonant signal is dominated by the structural anisotropy in the 
distribution of the F neighbors of the resonant Cu atom, and that the Cu-3d orbital ordering has 
only a minor influence on the spectra. Our LDA+U results, however, also indicate that a change 
in the magnetic structure has a small influence on the Jahn- Teller distortion, and hence on the 
resonant spectrum, in the conventional (room-temperature) crystallographic structure of KCUF3. 
This may indicate that the large change observed experimentally in the resonant signal near the 
Neel temperature is related to a low-temperature structural transformation in KCUF3. 

PACS numbers: 71.20.-b,71.70.Ch,78.70Ck 



I. INTRODUCTION 

The pseudocubic perovskite KCUF3 is a magnetic in- 
sulator that has attracted significant interest since the 
sixties, as a prototype material for orbital ordering, 
cooperative Jahn- Teller distortion, and low-dimensional 
magnetismi^ KCUF3 is also structurally related 
to high-T c superconducting cuprates and colossal- 
magnetoresistance manganites. In recent years, renewed 
attention has been given to the problem of orbital or- 
dering and to the study of the interactions between 
orbital, magnetic, and structural ordering in strongly 
correlated 3<i-transition-metal compounds, stimulated in 
large part by the interest in colossal-magnetoresistance 
manganites. 4 In the latter systems, the interplay be- 
tween orbital, charge, spin, and lattice degrees of free- 
dom is known to play an essential role in establishing 
the observed physical properties. Contrary to the case of 
charge, spin, and lattice ordering, however, which can be 
investigated by conventional x-ray, neutron, and electron 
diffraction techniques, direct experimental observation of 
orbital ordering has been — and remains in most cases — a 
challenge. 

Recently, resonant elastic x-ray scattering (RXS) tech- 
niques have been developed to probe orbital order- 
ing. RXS has been applied to a number of orbitally- 
ordered compounds, including various manganites^ and 
KCuFaiS Except for very recent measurements car- 
ried out on layered/doped manganites at the Mn £2,3 
edgggjSjiOjiiji^ mos ^ f the experiments performed to date 
have used excitations from the transition-metal Is core 



level, i.e., K-edge excitations. In the case of KCUF3 
and for related perovskite manganites such as LaMnC>3, 
the K edge excitation is the only possible core-level ex- 
citation for which resonant photons have a sufficiently 
large wavevector to access a Bragg vector of the tar- 
getted orbitally ordered structure. The K-edge scat- 
tering process in manganites and KCuF 3 , however, de- 
rives from virtual electric dipole excitations to empty 
transition-metal Ap states. The resulting RXS ampli- 
tude is thus not directly related to the 3d states that 
exhibit orbital ordering, but to the anisotropy induced 
in the surrounding 4p states. Although the latter is due, 
in principle, to the orbital ordering and/or associated 
Jahn- Teller distortion, this indirect dependence compli- 
cates the interpretation of K-edge RXS data. In fact, 
it has been first proposed that the measured anisotropy 
is determined by the 3d orbital polarization via onsite 
3d-Ap Coulomb interaction. 13 However, subsequent local- 
density-approximation (LDA) calculations and LDA plus 
onsite- Coulomb interaction (LDA+U) calculations, per- 
formed for a number of manganitesi 4 * 1 ^* 1 ^* 1 ^ and also for 
KCUF3 have provided rather strong evidence that the 
K-edge scattering amplitude in these systems is mostly 
sensitive to the Jahn- Teller distortion that accompanies 
the orbital ordering, rather than to the orbital ordering 
itself. 

In KCUF3, an interesting feature has been observed 
concerning the temperature dependence of the K-edge 
RXS intensity for orbital Bragg peaksi^ The measured 
scattering intensity shows a substantial and sudden in- 
crease (within 4-5 K), with decreasing temperature, when 
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approaching the Neel temperature. A somewhat similar 
phenomenon has been observed also in LaMnC>3 at the 
K edge£ and in Lao.5Sri. 5 Mn04 at the L2.3 edges^ for 
temperatures approaching the Neel temperature, but in 
these systems the changes in the orbital scattering inten- 
sity take place over a much wider temperature range, i.e., 
one to two orders of magnitude wider. These phenom- 
ena are therefore not necessarily related, but they are 
all very puzzling, because they take place in a tempera- 
ture region where both the orbital-order parameter and 
the Jahn- Teller distortion are expected to be saturated. 
In all of these systems, the origin of the changes in the 
scattering intensity is an unresolved issue. 

For KCuF 3 , two different interpretations have been 
proposed. On one hand, the change in RXS intensity has 
been ascribed to a change in the orbital order param- 
eter, and interpreted as indicative of a strong coupling 
between spin and orbital degrees of freedomi 7 ^ On the 
other hand, based on the results of LDA + U calculations 
showing a dominant influence of the Jahn- Teller distor- 
tion on the RXS spectra, the jump in intensity has been 
attributed to an increase of the Jahn- Teller distortion 
in the low-temperature magnetic phased These inter- 
pretations are also not necessarily mutually exclusive; a 
change in the magnetic structure may induce a change 
in the orbital-order parameter, which in turn may mod- 
ify the Jahn- Teller distortion. In principle, first-principle 
calculations could help discriminating between these var- 
ious possibilities by providing quantitative information 
on the influence of the magnetic structure on both, the 
Jahn- Teller distortion and the orbital-order parameter. 
As yet, however, to our knowledge these issues have not 
been addressed by such calculations for KCUF3. 

In the present work, we use state-of-the-art LDA+U 
calculations, based on ab initio pseudopotentials, to in- 
vestigate the electronic, structural, and magnetic prop- 
erties of KCUF3. Our main targets are the orbital or- 
dering and Jahn- Teller distortion, and we address the 
influence of the magnetic structure on these properties. 
We also evaluate the Cu K-edge RXS spectrum, and 
investigate the influence of the magnetic structure and 
Jahn- Teller distortion on the simulated spectrum. Unlike 
other density-functional-related approaches, the LDA+U 
method has not been extensively used to examine struc- 
tural properties. In this study we therefore also first 
demonstrate the ability of LDA + U pseudopotential cal- 
culations to provide a realistic description of the ground- 
state properties of KCUF3, when full structural optimiza- 
tion is performed. 

KCUF3 crystallizes in a perovskite structure which is 
tetragonally distorted by a collective Jahn- Teller distor- 
tion of the CuFg octahedra. Two different polytypes 
have been identified experimentally at room tempera- 
ture, which differ in their atomic-plane stacking along the 
c axis (Fig. QJ. 19 The Cu ions in KCUF3 are nominally in 
a Cu 2+ (d 9 ) configuration, with completely filled 3d-t 2g 
orbitals and one hole in one of the cubic degenerate 3d-e g 
orbitals (3d z 2_ x 2, 3d 3y 2_ r 2). The cubic degeneracy of the 
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FIG. 1: (Color online) Schematic views of the atomic struc- 
ture and of the ordering of the Cu-3d hole orbitals in the a- 
(left) and d- (right) polytypes of KC11F3. 



e g states is lifted in the paramagnetic phase, well above 
room temperature, by an orbital polarization and Jahn- 
Teller instability. The octahedra elongate along the a or b 
principal axes, with an alternate order on neighboring Cu 
A and B sites in the ab plane (see Fig.QJ. This distortion 
corresponds to an alternate distribution of 3d z 2_ x 2 and 
3d z 2_ y 2 hole orbitals on Cu A and B sitesi^ Along the c 
axis, the stacking of the planes may be antiferromagnetic- 
like (type a) or ferromagnetic- like (type d), giving rise to 
the two different polytypes. 

Below the Neel temperature (38 K for polytype a, 22 K 
for polytype d) , the spins arrange in an A-type antiferro- 
magnetic (AF) structure, with strong AF superexchangc 
interaction between adjacent Cu atoms along the c axis, 
and weak ferromagnetic (FM) superexchange coupling 
between Cu nearest neighbors within the ab plane^SS 
This magnetic structure is consistent with the predictions 
of the Goodenough-Kanamori-Anderson^ rules for the 
superexchange interactions associated with the orbital- 
ordered configurations displayed in Fig. ^ Such orbital 
configurations have been confirmed by previous LDA+U 
and Hartree-Fock calculations for KCuF3i 8 i ls i 21 i 22 i 23 Pre- 
vious calculations have also shown that the LDA ap- 
proach fails to correctly predict the stable orbitally- 
ordered AF insulating structure of KCUF3, while the 
LDA+U corrects this featured 

In this paper, we concentrate on the a-type structure 
of KCUF3, which has been the focus of the experimen- 
tal RXS studies. 7,8 The work is organized as follows: In 
Section [Bj we briefly present the calculation method. 
In Section IIIII we address the structural properties of 
the low-temperature A-AF structure; the corresponding 
orbital ordering and electronic properties are examined 
in Section IIVI In Section [V] we investigate the ener- 
getics and superexchange couplings of several different 
(meta)-stable spin structures of KCUF3. The influence 
of the magnetic structure on the orbital ordering and 
Jahn- Teller distortion is addressed in Section IVT1 In Sec- 
tion lVIII we present our results for the RXS spectra, and 
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we summarize our conclusions in Section IVIIII 



II. METHODOLOGY 

The calculations are performed within the LDA + U 
scheme, using the pseudopotential plane-wave method. 
We employ the rotational invariant LDA + U energy 
functionali 24 i 2 ^ 

E^^[p\p^{n^},{n^}} =£ LD Vy] + 

I,A,<7 

where E LBA [p^ , p^\ is the standard LDA energy func- 
tional, p a is the electron spin density with spin polar- 
ization a, U (J) is the onsite Coulomb (exchange) inter- 
action parameter for the localized electrons — the Cu-3<i 
electrons, and n^' CT are the occupation numbers of the lo- 
calized orbitals. I and A are indexes for atomic sites and 
orbitals, respectively. The occupation numbers n 1 ^ 17 are 
the eigenvalues of the density matrices: 

= E fU^WMAvti), (2) 

h.i 

where (p^f(r) are the rotational-degenerate atomic or- 
bitals of the localized-d-electron shell, m stands for the 
angular quantum number, V'fcil 7 *) arc the crystal Bloch 
states, and are the corresponding occupation num- 
bers. 

In our pseudopotential calculations, the crystal Bloch 
states, ip%i(r), are represented by Bloch pseudo wave- 
functions, and the matrix elements (¥?m <T |V'fe i)j m Eq. 0), 
are evaluated by projecting the crystal pseudo wavefunc- 
tions on atomic pseudo wavefunctionsj2£ We employ the 
Perdew-Zunger parametrization of the LDA exchange- 
correlation potentials 2 ^ For the effective onsite Coulomb 
interaction parameter, U e g = U — J, we use a value of 
6.6 eV. This value corresponds to the calculated values 
U = 7.5 eV and J = 0.9 eV obtained by constrained LDA 
computations for KCuF 3 in Ref. |^ 

We employ a Vanderbilt 28 ultrasoft pseudopotential 
for Cu and norm-conserving Troullier-Martins^ pseu- 
dopotentials for K and F. Some of the calculations have 
also been performed using a Troullier-Martins pseudopo- 
tential for CuiS 9 . for comparison. The pseudopotentials 
we use are scalar relativistic. The Vanderbilt pseu- 
dopotential has been generated in the non-magnetic Cu 
3d 9 5 As 5 Ap 1 atomic configuration, using for the core- 
cutoff radii: r s = 2.1, r p = 2.2, and r<j = 2.0 a.u.. 
The Troullier-Martins pseudopotentials have been gen- 
erated in the ground-state configuration of the non-spin- 
polarizcd atom, with the cutoff radii (in a.u.): r4 S = 3.50, 
rip = ?'3d = 3.75 for K and r 2s = r 2p = 1.40 for F. 
The Troullier-Martins pseudopotentials have been cast 
into the Kleinmann-Bylander'?- non-local form using the 



s-component as local part. For Cu and K, we have in- 
cluded the non-linear core correction^ to the exchange- 
correlation potential to account for the overlap between 
the valence and core charge; for fluorine this is not neces- 
sary, since we treat the F 2 s semicore electrons as valence 
electrons. 

All calculations are carried out in a 20-atom tetrago- 
nal unit cell. We use a kinetic- energy cutoff of 80 Ry 
for the plane-wave expansion of the electronic states. 
The integrations in reciprocal space are performed us- 
ing a (4,4,2) Monkhorst-Pack^ k-point grid for all struc- 
tures, except the non-magnetic structure of KCUF3 
which is metallic. For the latter structure we employ 
a (6,6,4) Monkhorst-Pack grid together with a Gaussian 
electronic-level broadening scheme,— with a full width at 
half maximum of 0.01 Ry. With the above parameters, 
the estimated convergence on the relative energies of the 
various magnetic structures considered in this work is 
~ 1 meV per formula unit. 

Concerning the pseudopotentials, we have performed 
some comparative studies using the Troullier-Martins 
and Vanderbilt pseudopotential for Cu. They indicate 
negligible differences, in general, for results concerning 
the electronic-structure properties. For the structural 
properties, the Vanderbilt pseudopotential yields a larger 
(~ 15 % larger) Jahn- Teller distortion, and a smaller 
equilibrium volume (2 % smaller) than the Troullier- 
Martins pseudopotential. This is related to a higher 
degree of localization within the core region of the Cu- 
3d pseudo-charge density obtained with the Vanderbilt 
compared to the Troullier-Martins pseudopotential. The 
structural trends, however, and in particular the depen- 
dence on the magnetic state, are the same with the two 
types of pseudopotentials. The effect of the K 3p and 3s 
semicore electrons has also been investigated. Potassium 
is known to have semicore states which are rather po- 
larizable. We have therefore performed test calculations 
treating the K 3p and 3s electrons as valence electrons. 
They show that the inclusion of the K semicore electrons 
has virtually no effect on the electronic properties and 
changes the equilibrium values of the structural parame- 
ters by less than one percent. 



III. STRUCTURAL PROPERTIES 

In Table [fl we show our results for the equilibrium 
lattice constant, a, tetragonal lattice ratio c/a, fluo- 
rine coordinate xp, and bulk modulus Bq of the low- 
temperature A-AF structure of KCUF3. The theoretical 
values of the structural parameters are compared to ex- 
perimental values at low temperature — when available — 
and at room temperature. The theoretical equilibrium 
lattice constant a is about 2 % smaller than the mea- 
sured low-temperature value. This is comparable to LDA 
results for related perovskite materials, such as the fer- 
roelectric oxides studied in Ref. 0i f° r which the LDA 
underestimation of the lattice constant is typically ~ 1 
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a (A) c/a x F Bp (GPa) 



LDA + U 


4.03 


0.955 0.2320 


75 


Expt a (10 K) 


4.126 


0.9486 




Expt 6 (298 K) 


4.141 


0.9476 0.2280 





TABLE I: Calculated values of the equilibrium lattice con- 
stant, a, tetragonal lattice ratio, c/a internal structural pa- 
rameter, xjj, and bulk modulus, B , of KCuF 3 . The experi- 
mental data are from Ref. [2f] (a) and from Ref. [3^ (b). 



to 2 %. 

Our LDA+U values for c/a and Xf agree well with 
experiment, i.e., to within 1 and 2%, respectively. The 
calculated relaxed bond lengths for the three inequiva- 
lent Cu-F bonds along the x, y, and z directions (at a 
Cu A site) are 1.87, 2.16 and 1.92 A, respectively, versus 
1.88, 2.25 and 1.96 A in the experiment^ The in-plane 
quadrupolar distortion AX = (0.25 — xf)2o, associated 
with the Jahn- Teller distortion of the CuF6 octahedra, 
is 3.6 % of a, compared to the experimental value of 
4.4 Hartree-Fock calculations 22 performed for the d- 
type structure yield a distortion AX /a of 2.6 %; the same 
Hartree-Fock calculations corrected a posteriori using a 
gradient-corrected correlation functional give a distortion 
of 3.4%, closer to our calculated value. Earlier LDA + 
U calculations, based on the full-potential linear-mufhn- 
tin-orbital method and performed using the experimental 
values of lattice constants, reported a smaller quadrupo- 
lar distortion of 2.5% for the cJ-type structured To our 
knowledge, no experimental value is available for the bulk 
modulus of KCUF3. Our value of 75 GPa, however, is in 
reasonable agreement with the -Bo value of 84 GPa de- 
rived from semiempirical pair-interaction-model calcula- 
tions for the elastic constants of KCUF3. 36 

The structural optimization needed to determine the 
parameter values in Table [3 was carried out consid- 
ering several different unit-cell volumes, and for each 
volume several different c/a ratios, and by fully relax- 
ing in each case the F internal coordinates. As a by- 
product, we have therefore access to the dependence of 
the quadrupolar distortion on tetragonal strain and also 
to the hydrostatic-pressure dependence of c/a and xf- In 
Fig. we display the calculated equilibrium quadrupolar 
distortion AX as a function of the c/ a ratio, for three dif- 
ferent volumes (including the experimental one) . Consis- 
tent with semiempirical model descriptions^ an increase 
in AX decreases c/a; the results in Fig. [^indicate that a 
10% increase in AX decrease c/a by ~ 5%. In Table ITT1 
we also reported the equilibrium values of Xf and c/a at 
different volumes. Decreasing the volume, decreases the 
quadrupolar distortion, and increases the c/a ratio. The 
latter trend is not inconsistent with the increase in the 
experimental c/a ratio observed with decreasing temper- 
ature (and volume) in Table U 

In the following sections, we report our results for the 
electronic properties, magnetic structure, and RXS spec- 
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0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 
c/a 

FIG. 2: Calculated quadrupolar distortion AX = (0.25 — 
XF)2a as a function of the cell-shape ratio c/a for different 
volumes of KCUF3. The arrows indicate the equilibrium val- 
ues under hydrostatic conditions. 



V (A 3 ) C/A X F 

67.03 0.937 0.2257 

63.68 0.951 0.2305 

60.39 0.961 0.2343 

57.24 0.968 0.2371 



TABLE II: Calculated equilibrium values of the tetrago- 
nal lattice ratio c/a and internal structural parameter if of 
KCUF3 at different volumes. 

tra of KCuF 3 evaluated at the reference experimental 
values of the structural parameters. We note that using 
the theoretical values of the structural parameters would 
not change significantly the results, i.e., produce negligi- 
ble changes in the orbital-order parameters, variations of 
~ 0.01/ib in the magnetic moments, and changes of the 
order of 0.1 eV in the relative energies of the electronic 
states. 



IV. ELECTRONIC STRUCTURE AND 
ORBITAL ORDERING 

In Fig. [31 we present the KCUF3 differential charge den- 
sity for the A-AF structure, i.e., the difference between 
the crystal electronic charge density and the superposi- 
tion of spherical atomic charge densities of neutral atoms. 



The differential density provides information on the elec- 
tronic charge rearrangements involved in the crystal for- 
mation, including charge transfers associated with bond 
formation and orbital ordering. The differential charge is 
displayed both in a CuF (001) [Fig. Eta)] and in a CuF 
(OfO) [Fig.Citb)] atomic plane of KCuF 3 . 

The differential density is clearly partitioned into F- 
related/Cu(K)-related electron accumulation/depletion 
regions, consistent with the ionic nature of KCUF3. The 
Cu-3d orbital ordering shows up quite strikingly in the 
maps, i.e., in the alternate sequence of depleted d z i_ x i- 
and d z 2_ y 2-like orbital densities on neighbouring Cu A 




FIG. 3: Difference between the crystal electronic charge den- 
sity and the superposition of spherical atomic charge densities 
of neutral atoms in a (001) (panel a) and a (010) (panel b) 
CuF atomic plane of KCUF3. The positions of the Cu (F) 
atoms are indicated by small (large) dark disks. Regions of 
electron accumulation (depletion) are shown in gray (white). 
Contour lines are separated by 0.03 |e|/bohr 3 ; negative elec- 
tron densities are truncated below -0.15 |e|/bohr 3 (near the 
Cu atoms). 
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FIG. 4: Electronic density of states of A-AF KCuF 3 . The 
zero of energy corresponds to the valence band maximum. 

and B sites in Figs. Efa) and (b). We notice a slight 
asymmetry in the d z 2_ x 2-X^&& (d z 2_ y 2-like) orbital den- 
sity related to the tetragonal distortion of the lattice, 
i.e., stronger depletions for the lobes along the x (respec- 
tively y) axis than for the lobes along the z axis. We 
also observe a local increase in the electronic density of 
some of the nominally "filled" Cu-3d orbitals — in partic- 
ular the d Sy 2_ r 2 (d 3x 2_ r 2) orbitals — suggesting a slight 
contraction of the 3d states of the Cu ion in the crys- 
tal with respect to those of the isolated atom. We note 
that the relaxation of the three inequivalent Cu-F bonds 
follows the trend of decreasing bond length with increas- 
ing Cu-3d localized-hole charge on the bond, as expected 
from electrostatic interaction with the F anions. 

The differential maps indicate a significant interaction 
(hybridization) between Cu-3g? e g and F-2p CT states and a 
related non- negligible derealization of the Cu hole on the 
nearest-neighbor F 2p a states. These features, together 
with the deviation from purely d z 2_ x 2 (d z 2_ y 2) hole or- 
bitals, suggest a somewhat incomplete z 2 — x 2 / z 2 — y 2 
orbital ordering. To be more quantitative in this state- 
ment, we have evaluated, from the density matrix, the to- 
tal number of 3d holes per Cu site, n^f e = 10 — J2a \ n A> 
and the expectation value of the number of holes in the 
d z 2_ x 2 orbital on an A site (or equivalently, d z 2_ y 2 or- 
bital on a B site), n^°'| = 2-^ CT n ^-x 2 ,z 2 -x 2 ■ We note 
that the occupation numbers we obtain from the density 
matrix for the ti g orbitals are identically 1.00 per spin, so 
that the total number of 3c? holes is identical to the num- 
ber of e g holes in our calculations. The calculated values 
are: n\f e = n h ° le = 0.63 and n A i hol t = 0.61. Hence, the 
degree of z 2 — x 2 / z 2 — y 2 orbital ordering, measured by 
- j °e j g s ig n i ncan ^iy reduced (by 39 %) with respect to 
its nominal value of one. Most of this reduction, however, 
can be accounted for by the decrease in the number of e g 
hole, n^° le , relative to one (37 %), and derives thus from 
the derealization of the 3d hole (due to hybridization) 
on neighboring F orbitals. The d z 2_ x 2 (d z 2_ y 2) character 
of the ordered orbitals on the A (B) sites, measured by 
the ratio n^i 1 ^. / n e° le : is instead nearly complete [i.e., 
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FIG. 5: Projected Cu-3d and F-2p atomic densities of states 
in KCUF3. Density of states projected onto: (a) the Cu 3d e 3 
(dashed line) and tig (solid line) orbitals, (b) the Fi, and (c) 
the F2 2p„ (dashed line) and 2p n (solid line) orbitals. In each 
panel, the upper half corresponds to the majority spin states 
and the lower half to the minority spin states. 



97% of the d z 2_ x 2 (d z 2_ y 2) type]. 

The calculated electronic density of states (DOS) of 
A-AF KCuF 3 is shown in Fig. 0] KCuF 3 is an insulator 
in our calculations with a gap of 1.54 eV. The valence 
DOS of KCUF3 exhibits two main structures: a low- 
energy structure (-7 to - 3.5 eV), which includes mostly 
Cu-3d states, with an admixture of F-2p states (mainly 
2p a states), and a high-energy structure (-3.5 to eV), 
which is dominated by F-2p states, and includes also a 
non-negligible contribution from Cu-3g? states. The cor- 
responding projected Cu-3g? and F-2p atomic DOS's are 
displayed in Fig. El^i The Cu-3d local atomic DOS has 
been split into contributions from the ti g and e g states 
[Fig. [3] (a)] and the F-2p local DOS's into contributions 
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FIG. 6: Density of states projected onto the Cu Ap x (solid 
line), Ap y (dashed line), and 4p z (dotted line) orbitals at a Cu 
A site. 



from the p a and p T orbitals on each of the two inequiva- 
lent F sites [Figs. [3] (b) and (c)]. The sharp lowest-energy 
conduction-band feature, in Fig. 0] clearly derives from 
the ordered Cu-e ff empty orbitals, that have the minority 
spin on each Cu site [see Fig. 0(a)]. Within the LDA+U 
scheme, the lowest-energy conduction-band feature may 
be viewed as the upper Hubbard band associated with the 
localized e g orbitals illustrated in Fig.^ the correspond- 
ing lower Hubbard band being the sharp lowest-energy 
valence-band feature in Fig. [5] (a). The existence of a 
significant hybridization between Cu-e ff and F-p a states, 
in Fig. is apparent from the energy degeneracies and 
analogous dispersions of the related features. Inspection 
of Fig.[S]also reveals a non-negligible interaction between 
Cu-t2 g and F-p n states. 

At higher energy, in Fig. 01 one finds a second 
conduction-band edge (at about 5 eV) associated with 
K-4s states, with also a contribution from Cu-4s states. 
The two sharp features located at ~ 9 and 11 eV, in 
Fig. 01 originate from K-3d states. We note that the 
Cu-4p states, which are the intermediate states in the 
Cu K-edge RXS process, contribute to the DOS over a 
large energy region that ranges from about 8 eV to 32 
eV (at least). The corresponding Cu 4p x , 4p y , and 4p z 
projected densities of states are represented in Fig. (for 
a Cu A site). The dominant features of the Cu 4p x and 
Ap v densities of states exhibit a splitting of about 6 eV 
produced by the F quadrupolar distortion AX. The in- 
ward/outward relaxation of the nearest-neighbor F an- 
ions along the x/y axis tends to shift the 4p x /Ap y states 
at higher/lower energy with respect to the 4p z states. 
It should be noted, concerning the high-energy unoccu- 
pied states occurring above 5 eV, that — unlike the lower- 
energy unoccupied Gu-d z 2_ x 2 states associated with the 
sharp feature — these states are totally insensitive to the 
effective Hubbard term, U e s, in our LDA+U calculations. 
These states are thus expected to suffer from the usual 
LDA underestimation of the gap (i.e., to be shifted to 



lower energy with respect quasiparticle energies obtained 
in GW calculations). 
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V. MAGNETIC STRUCTURES 

The Goodenough-Kanamori- Anderson (GKA) rules 
predict, 2 for the half-filled-d 2 2_ :E 2/(i 2 2_ a 2 orbital con- 
figurations displayed in Fig. ^ that the favorable Cu 
superexchange interactions should be antiferromagnetic 
and strong along the c axis, and weakly ferromagnetic 
within the ab plane. In order to probe the magnetic inter- 
actions in our system, and also to address the influence 
of the magnetic order on the electronic and structural 
properties (in the next section), we have considered sev- 
eral different possible (meta)-stable spin configurations 
for KCUF3, illustrated in Fig. [7| In addition to the A- 
AF structure, we have considered the AF G- and C-type 
configurations, the FM state, and also the non-magnetic 
(NM) configuration. 

In Table IIIII we give the relative energies of the vari- 
ous spin structures together with the corresponding val- 
ues of the ordered Cu spin moment m s — computed as 
the integrated absolute value of the magnetic moment 
per Cu atom. In agreement with experiment, and con- 
sistent with the GKA rules, the A-AF structure is the 
most stable configuration. The second-most stable con- 
figuration is the G-AF spin structure, where the AF order 
is preserved within the Cu chains along the c axis, but 
FM order has been replaced by AF order within the ab 
plane. The structures having, instead, FM order within 
the Cu c-chains (C-AF and FM in Fig. 0) are energet- 
ically much less competitive. For the latter structures, 
we note that AF spin order is more favorable than FM 
spin order within the ab plane (E[C-AF] < E[FM]). This 
is contrary to the situation observed for the structures 
having AF spin order within the chains, where FM order 
is preferred within the ab plane (E[A-AF] < E[G-AF]). 
With the exception of the E[C-AF] < E[FM] case, the 
energy ordering we find is consistent with the trend ex- 
pected from the GKA rules for half-filled orbitals. Con- 
cerning the E[C-AF] < E[FM] anomaly, we note that in 
our system the hole orbital occupation differs from the 
half- filled situation (see also next section), and for the 
weak superexchange interaction within the ab plane, such 
a deviation can change the sign of the superexchange cou- 
pling. 

The energy of the NM structure, in Table lllll is by 
far (0.2 — 0.3 eV) the largest of all. It should be noted, 
however, that a large fraction of the energy difference 
between the NM structure and the stable A-AF phase 
can be accounted for by the energy change of in isolated 
Cu 2+ ion when its configuration is switched from non- 
spin polarized to spin polarized [i.e., +0.24 eV in our 
LDA atomic calculations for a spin moment of 0.95 \xb 
(see Table III)]. For the paramagnetic phase thus, with 
disordered spins on the different Cu sites, a reasonable 
energy estimate would be: E[NM] -0.24 eV « 0.06 eV, 



A-AF G-AF C-AF FM 

FIG. 7: Magnetic structures considered in this work for 
KCuF 3 . 



which is intermediate between the energy of FM and A- 
AF phase. 

We note that the relative stability of the different 
structures, in Table IIIII is reflected in their bandgaps. 
We find that the bandgap changes from 1.54 eV, for the 
A-AF structure, to 1.47, 1.00, 0.72, and eV for the 
G-AF, C-AF, FM, and NM structures, respectively. 

In Table IIIII the ordered moment m s of a Cu with a 
given spin a increases, in the different spin structures, 
with increasing number of Cu nearest neighbours along 
the c axis having the same spin a. We also observe an 
increase in the moment, but significantly smaller, with in- 
creasing number of nearest neighbours having the same 
spin within the ab plane. These trends are consistent with 
the expected Cu-F-Cu spin-density overlaps associated 
with the orbital arrangements in Fig.^ The experimen- 
tal value of the Cu ordered spin moment measured by 
neutron diffraction in the A-AF structure at low temper- 
ature (4 K) is ~ 0.49/isr which is significantly smaller 
than our value of 0.91 fig. The small experimental m s 
in KCUF3, however, is generally attributed to zero-point 
fluctuations in the spin direction (beyond our calcula- 
tions), which lower the expectation value of the onsite 
spin moment. This effect is especially large in KCUF3, as 
a result of the quasi-one dimensional nature of its mag- 
netic structure. 

Based on the spin-polarized energies given in Table ITTT1 
one may derive some rough estimates for the Cu su- 
perexchange coupling constants along the c axis ( J c ) and 
within the ab plane (J a ). Using an Ising model — with 

interactions 2 J^^S^S^ between spins S = | on adja- 
cent A and B sites along the c axis (within the ab plane) — 
to map our calculated magnetic energies, we obtain for J c 
values of +70 and +48 meV from the energy differences 



E (meV) m s (fi B ) 



A-AF 0.91 

G-AF 10 0.90 

C-AF 58 0.98 

FM 70 1.03 

NM 297 



TABLE III: Calculated total energies and ordered Cu spin 
moments in the various magnetic structures of KCUF3 con- 
sidered in this work. Energies are given per formula units and 
measured relative to the energy of the A-AF structure. 



E[FM] - E[A-AF] and E[C-AF] - E[G-AF], respectively, 
and for J a values of —5 and +6 meV from E[G-AF] - 
E[A-AF] and E[FM] - E[C-AF], respectively. Clearly, the 
variation in the values derived for J c (J a ) indicates that 
the changes in the spin configurations that we consider 
are too large to be in the harmonic regime. These values, 
however, are nevertheless in order-of-magnitude agree- 
ment with the experimental values of J c — 17.5 meV and 
J a w -0.2 meV£2fi 



VI. INFLUENCE OF THE MAGNETIC 
STRUCTURE ON THE ORBITAL ORDERING 
AND JAHN-TELLER DISTORTION 



300 




A X (% of a) 



In Fig. |SJ we display the total energy of the various 
spin structures considered in the previous section as a 
function of the quadrupolar distortion AX. The calcu- 
lated equilibrium distortions range from 4.65 % (FM) to 
5.15% (NM), with an equilibrium value of 4.8% for the 
A-AF structure. The changes in the spin structure in- 
duce thus variations in the Jahn- Teller distortion which 
do not exceed 10%. This 10% value, obtained when 
the NM structure is considered, is a conservative upper 
bound for the estimated variation in the Jahn- Teller dis- 
tortion produced by the A-AF to paramagnetic trans- 
formation. Indeed, experimentally, strong antiferromag- 
netic spin correlations are known to persist above the 
Neel temperature within the Cu chains along the c axis. 
Locally, this would be more realistically described using 
configurations of the type G-AF and A-AF. For example, 
if one considers a structure composed of Cu AF chains 
along the c axis, with no interchain magnetic interac- 
tion, and uses an Ising model to describe the spin inter- 
actions, the energy of this structure is the average of the 
energies of the G-AF and A-AF configurations. Using 
such a description to model the structure above the Neel 
temperature, and with the energies of the G-AF and A- 
AF configurations given in Fig. |SJ the estimated change 
in the Jahn- Teller distortion produced by the magnetic 
transformation decreases to ~ 1 %. 

To compare the orbital ordering in the different spin 
structures, we have reported in Table llVl the calculated 
number of 3d z 2_ x 2 holes on the A site, n A j^° l %, together 
with the total number of 3d-e g holes per Cu site, n^ ole , in 
each structure. These values have been calculated using 
the same quadrupolar distortion for all structures, i.e., 
AX = 4.6 %. If one uses, instead, the calculated equilib- 
rium AX for each structure, the numbers, in Table ITV1 
change by less than 1 %. For all spin structures, the or- 
bital ordering is almost exclusively of the z 2 — x 2 / z 2 — y 2 
type; the d z 2_ x 2 /d z 2_ y 2 character, given by the ratio 

n ^ 2 '-a;2 / n e° le > ranges from 91 % for the NM structure to 
98 % for the G-AF structure. The major change we ob- 
serve is a strong decrease (~ 30 %) in the degree of orbital 
ordering, given by n^f^t , of the NM structure relative to 
those of the other structures. This decrease is associated 



FIG. 8: Total energy of the various magnetic structures of 
KCuF3 considered in this work as a function of the quadrupo- 
lar distortion AX. 



n A i hol l n\ ole 

A-AF 0.61 0.63 

G-AF 0.62 0.63 

C-AF 0.62 0.66 

FM 0.64 0.66 

NM 0.42 0.46 



TABLE IV: Calculated number of 3d-e g hole per Cu site, 
n eg , and expectation value of the number of 3d z 2_ x 2 hole on 
a Cu A site, n^' hole , in the various spin structures of KCuFs 

z2 — x 2 

considered in this work. 

with an increased delocalization of the 3d hole. Although 
there is an important change thus in the orbital ordering 
— and 3d hole occupation — induced by a change from 
the magnetic to the NM configuration, we will see in the 
next section that such a change has not a major impact 
on the RXS spectrum. 



VII. RXS SPECTRA 

Resonant elastic x-ray scattering is a second order pro- 
cess in which a core electron is virtually promoted to 
some intermediate states above the Fermi energy, and 
subsequently decays to the same core level. We con- 
centrate here on RXS near the Cu K-absorption edge, 
and examine the resonant scattering intensity for Bragg 
reflections that selectively probe the orbital order (or 
equivalently the Jahn- Teller structural order). The cor- 
responding Bragg vectors are G = (h,k,l), with h,k,l 
odd, in units (^, ^, -). We consider dipole transitions to 
Cu-4p band states, and neglect the effect of the core-hole 
potential on the intermediate states. 

The RXS intensity for a Bragg vector Q and incoming- 



(outgoing-) photon energy Hcu may be written as: 

2 



(3) 



where Rj are the positions of the Cu ions in the unit cell, 
e a {e' p) are the components of the polarization vector of 
the incoming (outgoing) photon, a(f3) = x,y, z, and 

^(M = 2^ - - F 4 P - — ^ ( 4 ) 

are the components of the resonant atomic scattering ten- 
sor for the jth Cu ion in the cell. In Eq. J3J|, ^/Jq is the 
Cu-ls wave function at site j, Eq is the corresponding 
core level energy, r(J) is the position operator measured 
relative to Rj, ip^ n (E^ p n ) are the empty Cu-4p band 
states (energies) with Bloch vector k and band index 
n, and T is the broadening corresponding to the inverse 
lifetime of the Cu ls-core-hole and 4p-electron excited 
states. 

For Bragg vectors associated with the orbital order, 
the structure factors, in Eq. |0, read: 
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where is the atomic scattering tensor of a Cu ion 

at site A (B), and the factor of two accounts for the two 
A (B) Cu ions in the unit cell. Because of the reflection 
symmetries a x , a v , and a z present at the Cu sites in 
KCUF3 , the only non- vanishing components of the atomic 
scattering tensors are: F^x , Fy^ B \ and F^} Z B \ and 
because of the 90-degree roto-translation symmetry that 
transforms the A and B site into each other, one has: 
F *Mv,v) = F vM^) and F t = F iz- The scattering 
intensity in Eq. J2J may therefore be written as: 



I(G,huj) oc \F£ x (Hu) - F A Jhw)\ \e x e' x - e y e' yl 



(6) 



Except for the special case where the polarization- 
dependent term on the right-hand side of Eq. © van- 
ishes, the scattering intensity is proportional thus to 
\F A x (hu>) — F y A y (huj)\ 2 . We also note that, within the 
extremely localized core region where the Cu-ls wave 
function does not vanish, the radial dependence of each 
band state il>u n , in Eq. is essentially the same as that 
of the atomic 4p state, except for a global scaling factor 
in the wave-function amplitude. The scattering factors, 
in Eq. ©, may therefore be written as: 



Fa,a{E) « / 
JE, 



deD 4 J>(e) 



E + Eq- e -iT/2' 



(7) 



where D^P(e) is the projected atomic 4p Q density of 
states and Ep is the Fermi energy. Hence, in our pseu- 
dopotential calculations we evaluate the orbital scatter- 



ing intensity, I or b(fiw) oc \F*r (fi/jj) 



F A 

v,v 



(M 



FIG. 9: Orbital RXS intensity (a) and absorption coeffi- 
cient (b), as a function of photon energy, near the Cu K edge 
in KCUF3. The calculated spectra of the A-AF (solid line) 
and NM (dashed line) structures are displayed. The effect 
of a 50% increase in the quadrupolar distortion, AX, in the 
A-AF structure is also shown (dashed-dotted line). The ex- 
perimental data are from Ref.Qf the measured RXS spectrum 
corresponds to a (3,3,1) orbital Bragg reflection and a a — n' 
polarization. 



function of photon energy, directly from the partial 4p 
densities of states. The core-level energy Eq (not com- 
puted here) is a value we adjust to align the main feature 
of the calculated and experimental RXS spectrum^ 

In Fig. Efa), we confront our calculated orbital RXS 
spectrum for the A-AF structure of KCUF3 to the exper- 
imental spectrum [measured at the (3,3,1) orbital Bragg 
reflection and for a a — it' polarization]. We also re- 
ported in this figure the theoretical RXS spectrum of the 
NM structure, for comparison. Both spectra have been 
computed using the experimental values of the structural 
parameters. A spectrum obtained for the A-AF structure 
with a 50% increase in AX is also shown in Fig. EJa). 
The corresponding calculated K-edge absorption spectra 
A(hu), evaluated as A(hw) oc Im[F A x (hw) + F A y (hjj) + 
F A z (huj)], are displayed in Fig. EHb), and also compared 
to experiment (fluorescence data) X All theoretical spec- 
tra have been evaluated using for T 0.5 eV, and we convo- 
luted these spectra with a Gaussian of full- width at half 
maximum of 1 eV. 
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The calculated RXS and absorption spectra for the A- 
AF structure are in reasonable agreement with the exper- 
imental data. We note, however, that the energy separa- 
tion between the two main features in the RXS spectrum 
is somewhat underestimated in our calculations, and the 
intensity of the high-energy feature near 9 keV is also not 
very well reproduced (both in the RXS and in the absorp- 
tion spectrum). The change in the magnetic configura- 
tion, from the A-AF structure to the NM structure, de- 
creases somewhat the RXS intensity. We note, however, 
that this change is very small compared to the factor-of- 
two decrease in intensity observed experimentally when 
the material transforms to the paramagnetic phase. The 
weak dependence of the RXS spectrum on the magnetic 
configuration is consistent with the results of previous all- 
electron calculations for the A-AF structure (calculated 
within the LDA+U) and the NM structure (calculated 
within the LDA in previous studies) 

To assess the influence of the orbital ordering on the 
RXS signal, we have also computed the RXS intensity of 
the A-AF structure without the quadrupolar distortion 
(AX ~ 0). Although the z 2 — x 2 / z 2 — y 2 orbital ordering 
is still large in this case: n 1 ° t = 0.64, n^ oZe = 0.68, the 
calculated RXS signal is about two orders of magnitude 
smaller than that obtained with the experimental Jahn- 
Teller distortion; the orbital ordering by itself thus has a 
negligible influence on the RXS intensity. We note that 
the fact that n 2 '_ 2 e (and n h ° le ) are very similar with 
and without the quadrupolar distortion indicates that, 
already at vanishing AX, the orbital ordering has essen- 
tially reached its saturation value, given the hybridiza- 
tion present in the system. 

Concerning the effect of the Jahn- Teller distortion, we 
observe that increasing AX by 50 % increases noticeably 
the splitting between the two mean features in the RXS 
spectrum of the A-AF structure, but does not change 
significantly the intensity of the dominant peak. De- 
creasing, instead, AX by the same amount has a larger 
influence on the RXS intensity^ but corresponds to a 
situation where the two features are essentially merged. 
A change in the Jahn- Teller distortion in the magnetic 
phase (not supported by our LDA+U calculations of the 
Jahn- Teller distortion) has been suggested as a possible 
cause of the drastic change observed experimentally in 
the RXS signal near the Neel temperature^ We observe 
here, however, that a factor of two change in the RXS 
intensity would require a drastic change in AX, in the 
conventional crystallographic structure of KCuF 3 . Such 
a change is bound to induce also a substantial change 
in the c/a ratio (see Section ILTl)l . which has not been 
observed experimentally. The Jahn- Teller origin of the 
jump in intensity appears thus somewhat unlikely. We 
therefore suggest that possibly a structural transforma- 
tion occurs at low temperature in KCUF3, which is re- 
flected in the RXS signal. We note that such a change in 
the crystal structure would not be inconsistent with some 
unexplained Raman- and nuclear magnetic resonance 40 
features observed experimentally at low temperature in 



this material. We also note that recently a new crystal 
structure (superstructure) has been proposed for KCUF3 
based on x-ray diffraction at room temperature! 41 " In our 
LDA+U calculations, however, his structure is found to 
relax to the conventional crystallographic structure of 
KCuF 3 . 



VIII. CONCLUSIONS 

We have studied by means of LDA+U pseudopoten- 
tial calculations, the structural, electronic, and magnetic 
properties of KCuF 3 , and investigated the Cu K-edge 
RXS spectrum for Bragg reflections associated with or- 
bital order. In our studies, we considered several different 
(meta)-stable spin structures for KCuF 3 in order to as- 
sess the influence of the magnetic structure on the orbital 
ordering, Jahn- Teller distortion, and RXS spectrum. 

For KCuF 3 , the LDA fails to correctly predict the sta- 
ble AF insulating structure. We have found here that 
LDA+U pseudopotential calculations provide an accu- 
rate general description of the properties of this system. 
The ground state is correctly predicted to be an A-type 
AF structure. The structural parameters obtained from 
the LDA+U calculations agree also well with experiment, 
with an accuracy comparable to that obtained by LDA 
calculations for other perovskite materials. The orbital 
ordering is predominantly of the z 2 — x 2 j z 2 — y 2 type, 
consistent with previous theoretical predictions for this 
system. We find, however, that the orbital-order param- 
eter is significantly reduced with respect to its nominal 
value due to Cu(3d)-F(2p) hybridization. We also find 
that, given this hybridization, the Cu-3d orbital order- 
ing in the A-AF phase is already saturated at vanishing 
Jahn- Teller quadrupolar distortion of the F neighbours. 

The RXS spectrum calculated for the A-AF structure 
agrees relatively well with the experimental spectrum. 
We find that the resonant signal is dominated by the 
Jahn- Teller distortion, with a minor influence of the or- 
bital ordering, in agreement with previous theoretical re- 
sults. Our LDA+U calculations, however, also indicate 
that a change in the magnetic structure has a small in- 
fluence on the Jahn- Teller distortion, and hence on the 
resonant spectrum, in the commonly accepted crystallo- 
graphic structure of KCuF 3 . We therefore suggest that 
the change observed experimentally in the RXS signal 
near the Neel temperature may be related to a low- 
temperature structural transformation in KCuF 3 . We 
also note that preliminary results obtained from LDA+U 
molecular-dynamic simulations for A-AF KCuF 3 indicate 
the existence of a superstructure, with tilted CuFg octa- 
hedra, which is slightly lower in energy than the conven- 
tional crystal structure of KCuF 3 . 
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